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1 Introduction 



The presence of a substrate favouring one of two coexisting bulk phases may 
cause this phase to form a film adsorbed at its surface. At a certain temper- 
ature T w below the bulk critical temperature T c the system may undergo a 
wetting transition at which the thickness of the film becomes macroscopic. 
Such wetting transitions have been extensively studied during the last years, 
especially for the planar substrates; for reviews see [1-3]. In addition, a lot 
of effort has been made towards understanding the modifications of system's 
wetting behaviour introduced by structured substrates, including chemically 
homogeneous nonplanar substrates [4-14]. In particular, for spherical and 
cylindrical geometries [15-25] the substrate's finite and positive radius R 
prevents the divergence of the adsorbed film's thickness so that wetting tran- 
sition point may be reached only in the limiting case R = oo. This is due 
to the fact that - unlike in the planar case - the area of the interface is not 
invariant under the change of its position - and increases with the increasing 
layer's thickness. The accompanying growth of the interfacial free energy 
prevents the wetting transition. 

To our knowledge not much attention has been paid to studying interfa- 
cial correlation functions in such systems, especially from the point of view of 
their temperature and substrate's radius depencence [16] . Our first aim is to 
systematically derive the effective interfacial Hamiltonian for adsorption on 
a disc including its radius and temperature dependence. This Hamiltonian 
is then used to evaluate the interfacial correlation function and the corre- 
sponding correlation angle parametrized by R and T. In particular, one is 
interested in seeing how the divergence of the correlation angle close to TV 
is prevented by finite values of R. 

In Section 2 we introduce the Landau-Ginzburg- Wilson theory for a 2- 
dimensional system consisting of two phases coexisting in presence of a 
disclike substrate. This theory serves to derive within the mean field ap- 
proach the effective interface Hamiltonian. We discuss the curvature-induced 
modifications of its form as compared to the planar case and propose a per- 
turbative procedure yielding approximate expression for the equilibrium in- 
terface position at a given thermodynamic state. 

Section 3 contains the discussion of the correlation function for an in- 
terface fluctuating around its mean circular shape. Within the Gaussian 
approximation the expressions for the correlation angle and the amplitude of 
the correlation function are derived. 

In Section 4 we apply the specific results obtained in Section 2 to evalu- 
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ate the correlation function. After introducing natural scaling variables we 
discuss the asymptotic behaviour of the correlation length and identify the 
different scaling regimes. In Section 5 we summarize our results. 



2 The Hamiltonian 



The standard method of deriving the effective Hamiltonians amounts to defin- 
ing a new set of variables describing the system's state on the coarse-grained 
scale and introducing constraints which specify a restiction on the allowed 
set of microstates compatible with fixed values of new variables. The ef- 
fective Hamiltonian is then extracted from the logarithm of statistical sum 
performed under this constraint. The mean field version of this procedure 
amounts to minimization of the microscopic Hamiltonian under the appro- 
priate constraint. 

Our derivation of the effective Hamiltonian for a 1-dimensional inter- 
face fluctuating in presence of a disclike substrate follows the lines first out- 
lined in Ref. [2Zj for the planar case. It starts from the description of a 
two-dimensional uniaxial ferromagnet based on the Landau-Ginzburg- Wilson 
Hamiltonian 
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where m(r) = m(r, <fi) is a scalar order parameter (magnetization), R - the 
substrate's radius, $(m) and $ 1 (m| r=fi ) - the bulk and surface free energy 
densities. The bulk free energy density at coexistence &o{ m ) exhibits two 
equally deep minima located at m a0 < and m^ > corresponding to two 
coexisting phases with opposite magnetizations. In our approach <&o(m) has 
the form of double parabola l2T ~ 



$ (m) 



—2~\m — m a0 ) 2 , 
^(m-m^o) 2 , 



for m < m* 
for m > m* 



(2.2) 



where K > 0, and p" 1 = £ 7 is the bulk correlation length in phase 7, 7 = 
a, (3. The parameter m* determines the value of the order parameter m at 
which the /?-phase goes over into the a-phase, i.e., it fixes the position of 
the a — (3 interface. The choice of m* value is arbitrary [21], pOj; in this 
paper it is set equal to zero [2Z1- The continuity of $0 at m = implies the 
relation Pama0 = — 1 between the model parameters. The double parabola 
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approximation to the general m 4 model makes further calculations tractable 
as it yields a linear equation for the constrained magnetization profile. The 
cost to be paid for this simplification is the necessity to limit the considered 
system's states to these far from bulk criticality (T <C Tc) (see |2E], [2*7]). 
For macroscopic disc radii this corresponds to p^R ^> 1 and £ 7 not too large 
as compared to the surface enhancement length ^ JH]- The opposite case 
was considered in Refs[18-19]. 

In the presence of external magnetic field h the free energy density $(m) = 
$>o(m) — hm may be rewritten as 

t / v j ^(rn-m a0h ) 2 + <£> a for m < 
$(m) = < K h (2.3) 
[ -2^(m - m m y + $^ for m > 0, 

where m j0h = m 70 + and $ 7 = —h (m^ + ) are the equilibrium bulk 
magnetization and free energy density of phase 7, respectively. Note that 
for h = one has $ 7 = 0. From now on we consider the thermodynamic 
states of the system corresponding to bulk coexistence; more precisely we set 
h = 0~ to stabilize the a phase in the bulk. 
The surface free energy density is reperesented by the expression 

$1(7%) = -/iirai - -gml, (2.4) 

where m\ = m(r = R, <ft), hi > is the surface magnetic field favouring the 
phase (3 at the substrate, and g is the surface enhancement parameter. In 
the present analysis we assume g < —Kpp which corresponds to a second 
order wetting transition of a reference planar system |27j . 

Minimization of Eq. (j2.1|) is performed under the crossing constraint m(r — 
l(4>), <p) = 0, where r = Z(0) > R specifies the position of the interface. This 
procedure yields the Helmholtz equation for the constrained equilibrium mag- 
netization fh 

KAm = -J (2.5) 
dm 

which has to be solved with the boundary conditions 



dr 



, (2.6) 



r=R dm l 

lim^oo m(r, / ) = m a0 , and lim r ^ r m(r, / ) = lim r ^+ m(r, l ) = 0. 
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Figure 1: A schematic illustration of an interface fluctuating around its equi- 
librium position I in presence of a disclike substrate with radius R. 



In the radially symmetric case, i.e., when I = Iq — const the constrained 
magnetization profile rho(r,lo) fulfils the zeroth order modified Bessel equa- 
tion. Its general solution has the following form 

m (r, Z ) - m~fO = C y K (p y r) + D 7 / (p 7 r), (2.7) 

where K and Iq are the 0-th order modified Bessel functions, while 7 = a 
corresponds to r > I, and 7 = (5 to r < I, respectively. The constants C 7 , 
are determined from the above boundary conditions: 

C a = -r^~ v D a = , (2.8) 



where 
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(2.10) 
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When the radially symmetric constrained magnetization is inserted into 
Eq. (j2.1j) one obtains - after subtracting the bulk terms - the following ex- 
pression for the effective interface Hamiltonian 



H[m ] = 27rV(R,l ) 
where the function V(R, Iq) can be split into three terms 

V(R,l ) = RXsp(R) + loZ a p{R,l ) + Ru(R,l ) 



(2.11) 



(2.12) 



The full expression for V(R, Iq) is presented in the Appendix. 

The quantity E^i?) represents the .R-dependent substrate-phase (3 sur- 
face tension, E ay g(i2, Iq) is the R and Zo-dependent interfacial stiffnes coef- 
ficient, and u(R,Iq) - the effective interface potential. The surface tension 
H s p(R) can be independently obtained by considering the disclike substrate 
immersed in the phase /3 and the explicit expression for it is presented in the 
Appendix. The radially symmetric equilibrium value of Iq will be denoted by 
I. It minimizes the function V(R,Iq), i.e., \i =i = 0. 

The next step is to find the constrained magnetization profile m(r, <f), I) 
in the case when small deviations from the radially symmetric case are 
present, i.e., m(r,<p,l) = fh (r, l ) + 5m(r, <fi,l) and l(4>) = + Sl((p). This 
can be done along the lines described in Ref. |27j. Within this approach 
5m(r,(f),l) = 9 ™°(* M 5l((f))(l + 0(5l 2 )). We are interested in the effective 
interface Hamiltonian evaluated relative to the radially symmetric equilib- 
rium configuration in which the interface has position r — I. Up to terms 
quadratic in d8l/d<p it takes the form 

.2 / \dd>) 



H Fl [5l] 



V(R,l) - V{R, I) 



(2.13) 



while the interfacial stiffness coefficient is identified as 

'<9mo(r, Z)\ 2 



^ a0 (R } l) = Kl 



dr- 

r 



dl 



(2.14) 



To proceed analytically we limit our considerations to the case of large 
substrate's radius p^R ^> 1, and to such choices of the thermodynamic state 
that Pj(l — R) 3> 1 and (I — R)/R <C 1. Then the dominant contributions to 
the above functions take the form 



T, a p(R,l) = a 



1 + °(p a l) + °(p3l\ 



+ 
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+ 2br 
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R<- , - S > 



+ O ((/ - R)e 



-3 Pf} {l-R) 



) (3.15) 
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lE a p(R,l) + Ru(R,l) = al 



+ ar 



(6r 2 + 



1 + 



PpRJ\ 



e -2 Pf3 {l-R) + q 




lR e ~P0(l-R) + 



(2.16) 



The coefficients a, a, 6 and c are expressed in terms of the L-G-W Hamilto- 
nian parametres: 



K 



a 



a = 2Kp l3 m l30l b = Kpp, 



c = KppQml^ s 20 = -Kp p (g 2 + 6Q - l)m 



where Q = — , and the parameter r 



hi+gmpp 



is interpreted as the 



distance from the planar substrate's wetting temperature (r ~ T ~ Tl-y ) . The 
positive values r > correspond to the wet state of the reference planar 
substrate. The expressions (|2.15j) and (|2.16|) reduce to those known for the 
planar case (see Ref . [2*7] ) upon neglecting terms proportional to powers of 
(I — R)/R as compared to one. 

As already remarked we are interested in system's states characterized 
by large values of R, and T close Tw- In the following we shall apply mean 
field approximation within which the equilibrium state minimizes V(R,lo)- 
Fluctuation effects missed in this kind of description are known to influence 
quantitatively the system's critical properties. For this reason our results 
will not give the correct indices in the planar limit [21]. They will also 
differ from those obtained in [T|3j, where the disc wetting phenomena were 
discussed within transfer matrix approach but no curvature modifications 
of the stiffness coefficient and the potential V were considered. We believe 
that the model discussed here correctly describes the curvature effects on the 
quantitative level. 

Consequently, the equilibrium thickness of the adsorbed layer X = I — R, 
fulfils the following equation 
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Its solution A can be determined approximately by considering the relevant 
case A <C R. The lowest order approximation A is obtained by neglecting in 
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Eq. (j2.17|) all terms proportional to 



^ - ppare-^ - 2 P(3 (br 2 + c)e- 2p ^ x ° = (2.18) 
R 



which leads to 



-or + JaV + 5^±£l 

For i? — > oo the above expression reduces to the well known solution for the 
planar case: pp\ (R = oo) = — log 2 {b^ T +c) • ^he approximation leading to 
Eq. (|2.19|) neglects all the curvature corrections to V(R,l) as compared to 
the planar case. To explore the curvature-induced effects we solve Eq. ()2.17|) 
replacing the expressions with In this way the following equation for 
Ai is obtained 

£ + + ^)- 1/2 e-^ - Pp ar(l + ^fV^ + (2.20) 

e ~ 2p e Xl = 



with solution 



-ar{\ + %)* + Jah*(\ + % ) + J^(&r» + c + c% 
PpK = - log — | - | Au 



4(6r 2 + c[l + f]) 

(2.21) 

By iterating this procedure one defines a sequence of approximations to A 
denoted A n , n = 0, 1, 2, . . .. Our analysis of the correlation function will be 
limited to n = 1. Solving numerically Eq. (j2.17j) and comparing the numerical 
result with the approximate solution Ai shows that the difference | (Ai — A)/A| 
does not exceed 10~ 6 (see 

It follows from Eq. (j2.21j) that the critical wetting state at which A diverges 
cannot be reached for finite R values. This confirms the well known fact of 
non-existence of wetting transitions for surfaces with positive curvature. 

3 Correlation function 



The interfacial correlation function is defined as 

<«(W«W>= r ns, -i3H Fl [si] • t 3 - 1 ) 



fVSh 
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For small gradients of the fluctuating interface and for small variations of its 
position relative to its equilibrium value 5/(0) = /(0) — Tone may approximate 
the Hamiltonian, Eq. (j2.13j) by the expression bilinear in 5/ 



H 



my 



(3.2) 



where V"{R,l) = ^ 

After applying the Fourier expansion to 5/(0) 

oo 

5/(0) = [An cos n<p + B n sin n<p} , 

n=0 

Eq. ()3.2)l can be rewritten as 
H Fl [5l] = \f:(Al + Bl)(n^^^ + V'\RJ))+lAlv'\RJ) 



(3.3) 



(3.4) 



71=1 



Inserting Eq. (j3.4|) into Eq. (j3.1|) . and integrating over A ,{A n ,B n } leads to 
the following expression for correlation function 



< 5/(0!)(5/(0 2 ) >= A{R,l) 



cosh 



i(l-|A0|/7T) 



sinh(^- 



(3.5) 



where A0 
A 



7T 



"V /(fl^i" * s ^ aB corre l a ti° n angle, and 
-r-77777T7r- - the correlation function's amplitude. The correlation func- 

/3V (R,l) 40 1 

tion depends on the absolute value of the relative angle A0 and is mini- 
mal at A0 = 7T. On the other hand, for the case A0 = one obtains 
W 2 =< 5/(0i)5/(0i) >= .4 coth(l/£0). In the planar limit R — > oo the 
quantity ^ corresponds to the parallel correlation length £|| divided by R. 

For further investigation of Eq. ()3.5|) one needs to specify I and V" (R, I) . 
In what follows we shall use the results of Section 2 by putting / = R + Ai. 
In this way one obtains 

V"(l,R) = 2vr[^ G r(i?+y)e- p ' 3Al +4^[(6r 2 + c) J R + cA 1 ]e- 2p ' 3Al ], (3.6) 

where the higher order terms have been neglected. The quantity H a p(R, /) is 
evaluated by substituting Ai into Eq. ()2.15|) . 
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4 Scaling regimes 



In this section we discuss the scaling behaviour of the correlation angle £^ and 
the amplitude A as functions of the variables r and R in the regime where 
|t| <C 1 and p^R ^> 1. We introduce the dimensionless scaling variables t 
and r 

ar a 2 t 2 ppR 

4c' T ~ 8ac 
such that Eq. ()2.19|) takes the following form 



t 



(4.1) 



P/3A0 = - log 



-t+A t 2 



t 2 



r J 



(4.2) 



Note that the variable r depending both on the magnitude of r and 1/ppR 
may cover arbitrary region of the positive semiaxis, whilst p 1 R 3> 1 implies 
r/t 2 > I. 



The leading asymptotic behaviour of the correlation angle £^ and the 
amplitude A is determined by Ao- In this case one has to distinguish four 
regimes in the section of the R — r plane such that ppR ^> 1 and |r| <C 
1 depending whether T > Tw, T < Tw, r 3> 1, and r < 1. However, 
introducing Ai instead of Ao, as we do, modifies the correlation angle ^ 
and the amplitude A on the level of corrections to the leading asymptotic 
behaviour. It turns out that from this point of view eight subsets of the 
asymptotic regime ppR 3> 1 and \t\ -C 1 must be distinguished. They are 
listed below together with the corresponding expressions for Ai, and A: 

1. T <T w ,r^> 1, |rt 2 log(-2t)| < 1 

pp\ x = -log(-2t)(l + i ^t y + H.O.T.) 
U = -&*-* l *yfil{l-\\ + B.O.T) (4.3) 
(3ap P A=^^f c \(l-l\ + H.O.T). 

2. T < T w , r > 1, |rt 2 log(-2t)| > 1 



r p p x, = -io g (-2t)(i + ^4^) + H -°- T -) 

\ ^ = -^ 3 / 2 yfi(l + H 2 log(-2t) + H.O.T.) (4.4) 
( Pop p A=^y^f c \(l + kt 2 \og(-2t) + £O.T.), 
where k = —4-. 
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3. T < T w , r <C 1, r 3 / 2 > |t 2 log ^| 

f PA = -log^(l + v^m^) + ^-O.T.) 

U = -K 3/2 yf - ^ + ^-O.T.) (4.5) 
[ (3ap,A = -\^f c ^[l - \V~r + H.O.T). 

4. T<TV, r«l,r 3 / 2 « |t 2 log ^| 

r P/J Ax = -log^(l + ff + TT.O.T.) 

^ = -K 3/2 yf ^(l + ^log-| + TJ.O.T.) (4.6) 
[ /fop^ = -Iv^Ix (l + log ^ + H.O.T). 

5. T>TV, r>l, |t 2 log£|<l 

f PA = -log| : (l- lfl ^ y + //.O.r.) 

U = ^- 3/2 ^{l + ^>zi + H-O.T.) (4.7) 

6. T>T W , r>l, |t 2 log£|>l 

f P/J Ai = -k>g£(l-2f£ + ff.O.T.) 

^ = >- 3/2 V?^(l + lrflo g ^ + ^O.T.) (4.8) 
( /^P/^ = f v^T (l + f^log± + H.O.T.). 

7. T > TV, r < 1, r 3 / 2 > |t 2 log 

r p / ,A 1 = -tog^(i- E5 ^ + ^.o.r.) 

^ = K 3/2 V^( 1 + ^ + #.O.T.) (4.9) 

[ ^^ = 1^^(1 + 5^+^-0.7.). 

8. T > TV, r < 1, r 3 / 2 < |t 2 log 

r p^-log^l + ff + lf.O.r.) 

l e^R^/f^^ + ^log^ + TT.O.T.) (4.10) 
1 ^P^=lxAlf (l + ^log^ + TT.O.T.). 
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In all the specified regimes of parameters R and r, the quantity Ai di- 
verges logarithmically either as function of r or as function of R. In the 
above regimes 1 and 2 the leading term is only temperature dependent, i.e., 
ppXi ~ — In |r| whilst for the cases 3, 4, 7 and 8 it depends exclusively on 
the radius, i.e., p@ Ai ~ lnfjogi?). The leading contribution to the correlation 
angle is proportional to jj- for regimes 1 and 2, which for the corresponding 
correlation length ^\\ = R^^ gives £|| ~ A. On the other hand, for cases 3-8 

we obtain £^ ~ i? -1 / 2 , which means that the dominant term is temperature- 
independent and converges to in the limit R — > oo. In these cases one 
has ^|| ~ R 1 / 2 . For the amplitude A we obtain A ~ r _1 in cases 1, 2, and 
A ~ R 1 / 2 in regimes 3-8. Note that in regimes 3-8 the interfacial width 
W becomes proportional to R 1 / 4 . We also observe that the dependence of 
the leading term to the derived expressions for £^ and A on the free a — j3 
interfacial stiffness coefficient a is not the same for all the above specified 
cases. In particular, at fixed R and r in regions 1 and 2 the quantity £^ is 
proportional to \fa while in regimes 3-8 the leading contribution to does 
not depend on a. On the other hand, as far as the dominant contribution is 
concerned, in regions 1 and 2 the amplitude A ~ 1/ \fo while in all the other 
cases - A ~ 1/cr. 

Among the derived corrections one encounters terms which are nonana- 
lytic in the zero curvature limit when the temperature approaches the wetting 
temperature of the reference planar system. (For discussion of non-analytic 
contributions to other thermodynamic quantities in wetting systems, see eg. 
Refs.|2II], [21] )• These terms contain a logarithm of R or r multiplied by a 
power of 1/R or r, which is characteristic for the considered potential V(R, I) , 
appropriate for short-range interactions. In our treatment of the correlation 
function the presence of the non-analytic terms is a consequence of the struc- 
ture of the potential V(R, I) and the stiffness coefficient S a( g(i?, I) derived in 
Section 2. The form of corrections to the leading asymptotic behavior of ^ 
and A follows from the curvature-dependent potential V(R,l) in regimes 1, 
3, 5, 6, and 7. In regimes 2, 4, and 8 their form follows exclusively from the 
necessity of incorporating the position dependence of the stiffness coefficient 
T^ a p(R,l) into the analysis. To consider still the next to leading order cor- 
rections and A one should use - within the present approach - a better 
approximation to A (say, A 2 instead of Ai). However, additional restrictions 
must then be imposed on the parametres r and R to determine the charac- 
ter of the next correction; accounting for it causes necessity to increase the 
number of regimes of parametres r and R that must be distinguished. 
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5 Summary 



Starting from the Landau- Ginzburg- Wilson theory we derived within the 
mean field approximation the effective Hamiltonian for the fluctuations of 
an interface separating two coexisting phases, one of them forming a layer 
adsorbed on a disclike substrate. The Hamiltonian is parametrized by the 
temperature T and the disc radius R, i.e., both the coefficient of the inter- 
facial stiffness and the interface potential V appearing in the effective 
Hamiltonian depend on these two parameters, and on the interface position 
I. We applied the Gaussian approximation and found the expression for the 
height-height correlation function, in particular the correlation angle £^ and 
the amplitude A. Considering the case of large disc radii and temperatures 
close to the wetting temperature Tw of a reference planar system and keep- 
ing track of the position dependence of the stiffness coefficient £ a/ 3 and the 
curvature corrections to the effective interaction potential V we identified 
the scaling regimes in the space of parametres R and T — Tw in which the 
dominant contributions to £^ and A were determined. In addition we have 
derived the expression for the .R-dependent surface tension coefficient S S/ g 
characterizing the substrate-phase j3 interface. This surface tension coeffi- 
cient can be represented as the power series in 1/R, and the corresponding 
Tolman's length [35], [31)] turns out to be proportional to (T — T w ) 2 . 

In the limit of large disc radii our results for the equilibrium adsorbate 
layer thickness A and the correlation length ^\\ ~ R^ differ from those ob- 
tained by Gil and Mikheev (G-M) in Ref.^HJ via the transfer matrix ap- 
proach. These authors' analysis referred to a temperature independent step- 
like potential Vs(X)/R = for A > 0, Vs(\)/R = oo for A < which captures 
the relevant features of complete wetting phenomena in the planar limit. For 
that case they predict I ~ R 1 ^ 3 and £\\ ~ R 2 ^ 3 . The discrepancy between our 
and G-M result is due to fluctuation effects missed in our description. For 
the same reason we fail to predict the correct values of the critical indices in 
the planar 2D limit [31 j . Performing a transfer matrix calculation taking into 
account the effect of curvature on the structure of the effective Hamiltonian 
is a demanding task. We believe that the MFT approach applied here gener- 
ally accounts for the observable phenomena, though obviously the system's 
behaviour is further modified when fluctuations are incorporated into the 
description. It is worth noting that the phenomenological treatment based 
on balancing the curvature-induced Laplace pressure and the disjoining pres- 
sure applied by G-M to the case of long range infractions yields in the case 
of short range forces the results consistent with ours. 

From this perspective it would be very interesting to analyze the interfa- 
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cial correlation function in system with short range forces by properly taking 
into account all the relevant fluctuations. It has been shown [23 EB] that such 
planar systems may undergo either first or second order wetting depending 
on the system parameters. However, this alternative applies to systems with 
dimension larger than d c = 2.41, [H3], which is not the case considered in this 
paper. On the other hand, whether the interfacial correlations in the case 
of cylinder- or spherelike substrates contain some indications of the possible 
change of the order of the wetting transion in the corresponding planar case 
remains an open question. 
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6 Appendix 



In this Appendix we quote the complete expressions for the effective interface 
Hamiltonian V(R, I) in the radially symmetric case and for the substrate- 
phase P surface tension S Sj g(i?) as derived from the L-G-W theory. These 
expressions have the following form: 



V(R,l) 



K 



Cffi l 2 K (ppl)(K ( Pf3 l) - K 2 (p p l)) - R 2 K ( Pp R)x 



x 



K (ppR)-K 2 (p p R))\ +-D}pl [l 2 I (ppl) [Io(ppl)-l2(ppl)) -R 2 I (ppR)x 
x (h{p p R) - hippR))] + ZCpDp [K {u)I {u) - K^hiu^udu 

V /J JppR 



+ 



+2nR 



and 



1 



- hmpoh - - gm poh - ( h+gmpoh J CpK (ppR) - ( h+gmpm ) Dph (ppR) + 



l ~gClKl(p p R) - ^gDplippR) - gCpDpK Q (ppR)h{p p R) 



+ 



K 2 (pJ) 1 
lK (pJ) 



(A.1) 



E sP {R) = -h m m - \gm 2 mh + 
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-2 9 



K /Fi(R)\ 2 r 
- pbRK (ppR) [K 2 {p p R) - K {p p R) 

(^|f|) 2 ^oM + (hi + 9 m m )^K ( Pp R) , (A.2) 



where the functions Fi(R) and i^C-R) are defined in the main text, see 
Eq. (j2.10|) . The integral J^ l R [K (u)I (u) — Ki(u)Ii(u)]du in the above expres- 
sion for V(R, I) can be performed analytically for the case ppR 3> 1, ppl 1 
by applying the asymptotic expansion of the Bessel functions. The asymp- 
totic expansion of this formula is quoted in the main text, see Eqs fl2.15|) . 

(f2~m 

The expression for the surface tension S S/ g(i?) can be derived by analyzing 
the the L-G-W theory for a disc immersed in phase (3. One notes that the 
function S S/3 (_R) can be expanded in powers of 1/ppR and the first two terms 
of this expansion have the following form: 

E S/3 (R) = r 2 - hi (mpo + r) - ^{mp + rf + 

^^- + H.O.T. . (A.3) 

Thus we observe that the quantity identified as the Tolman's length depends 
on the temperature and vanishes upon approching the planar substrate wet- 
ting temperature TV proportionally to r 2 . This behavior of the Tolman's 
length reflects the general property of Yi a p(R) which can be read off from the 
above formula, namely that at T = T\y the surface tension ceases depending 
on R. 



Finally, within the L-G-W theory the expression S a/ g(i?, /) for the R and / 
dependent interfacial tension is identified as T, af3 (R, I) = Kl dr^ 9 ™ 1 ^ ) 2 . 
The integration can be performed in the asymptotic regime p 1 R 3> 1 and 
the result is quotes in the main text, see Eq. (j2.15j) . 



15 



References 

[1] S. Dietrich in Phase Transitions and Critical Phenomena edited by C. 
Domb and J. Lebowitz (Academic Press, New York, 1988), Vol 12. 

[2] M. Schick in Liquids at Interfaces, 1988, Les Houches Lectures , edited 
by J. Charvolin, J. F. Joanny and J. Zinn-Justin (North Holland, Ams- 
terdam, 1990) 

[3] P. G. de Gennes, Rev. Mod. Phys. 57 (1985) 827. 

[4] R. R. Netz and D. Andelman, Phys. Rev. E 55 (1997) 687. 

[5] H. T. Dobbs and J. M. Yeomans, Mol. Phys. 80 (1993) 877. 

[6] W. R. Osborn and J. M. Yeomans, Phys. Rev. E 51 (1995) 2053. 

[7] K. Rejmer and M. Napiorkowski, Z. Phys. B 102 (1997) 101. 

[8] P. S. Swain and A. O. Parry, Eur. Phys. J. 4 (1998) 459. 

[9] C. Rascon, A. O. Parry and A. Sartori, Phys. Rev. E 57 (1999) 5697. 

[10] K. Rejmer, S. Dietrich and M. Napiorkowski, Phys. Rev. E 60 (1999) 
4027. 

[11] K. Rejmer and M. Napiorkowski, Phys.Rev. E 62 (2000) 588. 

[12] G. P. Kubalski and M. Napiorkowski, J. Phys.: Condens. Matter 12 
(2000) 9221. 

[13] A. O. Parry, M. J. Greenall and A. J. Wood, J. Phys.: Condens. Matter 
14(2001) 1169. 

[14] G. P. Kubalski, M. Napiorkowski and K. Rejmer, J. Phys.: Condens. 
Matter 13 (2001) 4727. 

[15] T. Bieker and S. Dietrich, Physica A 252 (1998) 85. 

[16] T. Gil and L. V. Mikheev, Phys. Rev. E 52 (1995) 772. 

[17] M. P. Gelfand and R. Lipowsky, Phys. Rev. B 36 (1987) 8725. 

[18] J. O. Indekeu, P. J. Upton and J. M. Yeomans, Phys. Rev. Lett. 61 
(1988) 2221. 



16 



[19] P. J. Upton, J. O. Indekeu and J. M. Yeomans, Phys. Rev. B 40 (1989) 
666. 

[20] R. Holyst and A. Poniewierski, Phys. Rev. B 36 (1987) 5628. 

[21] R. Evans, R. Roth and P. Bryk, to be published. 

[22] H. D. Wagner, J. Appl. Phys 67 (1990) 1352. 

[23] P. Levinson, J. Jouffroy and F. Brochard, J. Phys. Lett. 46 (1985) L-21. 

[24] I. Hadjiagapiou, J. Chem. Phys. 105 (1996) 2927. 

[25] F. Brochard, J. Chem. Phys. 84 (1986) 4664. 

[26] R. Lipowsky, Z. Phys. B 55 (1984) 345. 

[27] A. J. Jin and M. E. Fisher, Phys. Rev. B 47 (1993) 7365. 

[28] K. Rejmer and M. Napiorkowski, Z. Phys. B 97 (1995) 293. 

[29] G. Backx and J. O. Indekeu, Physica A 218 (1995) 69. 

[30] A. O. Parry and C. J. Boulter, J. Phys. Condens. Matter 6 (1994) 7199. 

[31] D. B. Abraham in Phase Transitions and Critical Phenomena edited by 
C. Domb and J. Lebowitz (Academic Press, London 1986), Vol 10. 

[32] The numerical check of the proposed approximation's accuracy was car- 
ried out for f = 2 = f = 1, ppR G [10 2 , 10 4 ] and r G [-0.5, 0.5]. 

[33] A. J. Jin and M. E. Fisher, Phys. Rev. B 48 (1993) 2642. 

[34] C. J. Boulter, Phys. Rev. E 57 (1998) 2062. 

[35] R.C. Tolman, J. Chem. Phys. 17 (1949) 333. 

[36] J.R. Henderson, Mol. Phys. 50 (1983) 741. 



17 



